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1. INTRODUCTION 

An electronic nose (e-nose) is a sensing instrument comprising a set of different gas sensors that 
react when exposed to a wide range of chemical particles. According to the outputs obtained from the 
sensors, many conclusions about air quality and toxicity can be obtained. They have shown great promise and 
utility in three domains: food control, disease diagnosis, and environmental monitoring. In food control, e- 
noses have been devoted toward assuring quality and safety for consumers. Some features already addressed 
in this domain are food freshness, ageing, contamination during processing, shelf life, and authenticity 
confirmation. Viejo et al. [1] implemented an e-nose to assess the aroma profiles in beer and automate the 
industrial quality inspection process. Similarly, Radi et al. [2] developed an e-nose for classifying odors from 
synthetic flavors such as grapes, strawberry, mango, and orange. To alert about rancidity, Xu et al. [3] 
introduced an e-nose monitoring the changes of pecans during storage. Timsorn and Wongchoosuk [4] 
explored the use of an e-nose device to identify odors from formalin contamination in seafood. To prevent 
products’ adulteration, Swigto and Chmielewski [5] proposed an e-nose to assist in the authenticity testing of 
products such as meat, honey, milk, and plant oils. To discourage meat dealers from committing food fraud, 
Laga and Sarno [6] presented an e-nose discriminating pork from beef. Wang et al. [7] deployed an e-nose 
inside a domestic refrigerator to assess the food freshness level of fruits, vegetables, and meat. 

The ability of humans to detect diseases with smell has played a significant role in clinical 
diagnosis. E-noses can ease the detection of volatile organic compounds (VOC) exhibiting bacterial 
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pathogens and have the potential to become a valuable disease diagnosis tool for humans, animals, and 
plants. Sanchez et al. [8] monitored exhaled human breath with an e-nose to help diagnose and monitor 
certain digestive and respiratory diseases. Siyang et al. [9] analyzed the odor of human urine samples with an 
e-nose to detect diseases such as diabetes. In response to the current COVID-19 pandemic, airbus is 
deploying e-nose devices in its aircrafts to detect the virus [10]. In veterinary medicine, Jia et al. [11] 
explored an e-nose device for detecting three types of wound infections in rats. To detect bovis-infected 
cattle in farms, Peled et al. [12] reported the use of an e-nose to analyze VOC in breath samples. In botanics, 
Baietto et al. [13] explored the use of an e-nose to sense the bole-rot fungi that affects trees. Wilson [14] 
developed and tested an e-nose for the rapid identification of insecticide residues in crops. 

Nowadays, e-noses in environmental monitoring have found application in four fields: i) air quality, 
ii) water quality, iii) process control, and iv) odor control systems [15]. Some representative work include 
Wongchoosuk’s WiFi e-nose sensing and quantifying indoor air contaminants even in very low 
concentrations [16]. The e-nose prototype proposed by Mishra devoted to identify poison gases emanated 
from waste [17]. The work of Baby reporting on the use of an e-nose to sense contaminating residues and 
pesticides in water [18]. The system in [19] using different gas sensors, monitors gas concentration and 
temperature in a biogas reactor. The smart system proposed by [20] was capable of sensing ammonium 
nitrate that could lead to fire and explosion in storage warehouses. Applications of e-noses outside these three 
main domains also involve explosive detection [21] and the space industry [22]. 

Within the context of air quality assessment, this paper presents the characterization and 
implementation of a novel, affordable e-nose device. The prototype consists of a set of six metal oxide 
semiconductor (MOS) gas sensors capable of detecting, i) combustible gas, ii) alcohol, iii) methane, 
iv) carbon monoxide, v) hydrogen, and vi) smoke, which might represent a threat for the human health. The 
implemented device has small dimensions allowing for seamless integration into mobile robots, indoor 
facilities, urban transport, risky environments (such as mines and sanitary landfills), or sensor networks 
monitoring broad surface areas. 

Similar e-nose prototypes exploring MOS sensors can be found in the literature [1], [16], [23]-[25]; 
sensor arrays range from four to ten, depending on the application. Prototypes rely on MOS technology 
because it offers small-size and robust sensors, quite good sensitivity, simple signal processing, commercial 
availability, and low cost. The main difference across devices is the method for data processing. Approaches 
such as artificial neural networks (ANN), principal component analysis (PCA), Fourier transform, and 
wavelets have been explored with satisfactory results for further developing predictive models. 

In this paper, the sensors’ signals were spectrally analyzed using the filter diagonalization method 
(FDM). To our knowledge, this work is the first one that reports on the use of FDM for e-nose signal 
processing. The remainder of this work paper is structured as follows: section 2 introduces the e-nose device, 
the MOS sensors used, the main system components, the FDM, and its implementation in the prototype. In 
section 3, the experiments that have been carried out are described and the results are shown. Finally, the 
concluding remarks and the future work perspectives are given in section 4. 


2. RESEARCH METHOD 

An operational prototype was implemented with low-cost materials, and it is the first approach for 
the research work. The proposed device involves six gas sensors, each one focusing on a specific gas. They 
capture the components present in the air for later off-line analysis. The device comprises an electronic 
module for data acquisition and a software to visualize the sensors’ behavior. To ensure an adequate gas 
concentration around the e-nose, two items have also been considered: an air pump conveying the samples 
and a hermetic box enclosing the device. 


2.1. Experimental apparatus 

A 3D printed plastic base was designed to host the set of sensors. It is small (6.5x5x7 cm) and 
displays six perforations around its side-faces for the proper installation of the MOS sensors as shown in 
Figure 1(a). All electronic boards and electrical connections are placed inside the plastic base allowing a 
clean and safe device handling as shown in Figure 1(b). 

The prototype employs six different types of MOS sensors (Hanwei Electronics Co., Ltd., Henan, 
China), which change their electric resistance when its sensing material comes in contact with the gas. 
Figure 2 illustrates the MOS sensors external and internal structures. The external structure as shown in 
Figure 2(a) involves a mesh-like enclosure that protects the sensing material and filters out the suspended 
particles on the air, so that only gaseous elements access into the chamber. The clamping ring secures the 
mesh and serves as the base for the sensing material and the electrical connections. 

Internally, the sensor is composed of a tin dioxide (SnOz2) layer, which is the sensing material that 
reacts to the input gas. In clean air, no electric current flows through this layer, but when gas is detected, 
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electrons are released, allowing current to flow. The current output signal is conveyed via platinum (Pt) wires 
as shown in Figure 2(b). Using a simple voltage-divider configuration, the gas concentrations can be derived. 
Each of the six MOS sensors comprised in the prototype is sensitive to a specific gas providing valuable data 
upon the analysis of their responses. We previously reported a complete experimental characterization of the 
MOS sensors [26]. Characteristics such as sensitivity, behavior to temperature, and step response were 
examined. Figure 3 presents, as an example, the behavior obtained from the alcohol (ethanol) gas sensor (S2). 
Figure 3(a) shows its resistance ratio Rs/Ro; here, Rs represents the sensor resistance to the target gas 
(C2HsOH) given a specific concentration, while Ro represents the sensor resistance in clean ambient air. Note 
that the resistance ratio Rs/Ro decreases as the concentration of the gas increases. Figure 3(b) shows how the 
sensitivity curve of the sensor is affected by temperature. Here, Rso represents the resistance of the sensor in 
125 ppm (parts-per-million) alcohol at 20°C. Note that the resistance ratio Rs/Rso decreases as 
environmental temperature increases. Figure 3(c) shows the sensor response to a step input of 125 ppm 
alcohol gas. Note that S2 delivers 3.5 V in the presence of the gas and 0.5 V (the baseline) in the absence of gas. 
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Figure 1. The plastic base structure: (a) 3D design and (b) actual prototype 
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Figure 2. MOS sensors: (a) external and (b) internal structures 


Table 1 lists the MOS sensors used, their model number, and their target detection gases. Note that 
the sensors are sensitive to a primary gas but also react to secondary gases, resulting in an overlap of their 
responses as shown in Figure 3(a) and Table 1. This overlap is advantageous since it provides more 
information about the chemical compounds contained in the air sample, thus facilitating the signal 
classification stage. 

The prototype also comprises an electronic module which encompasses a 32-bit microcontroller 
(Texas instruments Stellaris LM4F120) responsible for the sensors’ signal acquisition, universal serial bus 
(USB) interfacing with a computer, and the external air pump control. Figure 4(a) shows the block-level 
diagram of the electronic module while Figure 4(b) shows the electrical diagram for the air pump control. To 
synchronize the start and end of signal collection and transfer the acquired signals for processing and 
visualization, a simple communication protocol was built between the electronic module and the computer. 

Figure 5 shows the activity diagram of the system operation. The computer initiates the process 
which lasts 180 s. During this time interval, the sensor array collects data from its vicinity. At t=20 s, the 
pump starts introducing air into the hermetic box. At t=50 s, the pump is set off; the sensors reset (i.e., signals 
start going to their baseline). The electronic module captures this behavior as well. 
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Figure 6 illustrates all six sensors’ responses detailing the different stages. Note that the response of 
S2 (alcohol gas sensor) is the slowest one both in its rising and recovery times. This result is due to the low 
volatility that alcohol exhibits compared to the other substances explored. 
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Figure 3. Sensor characterization example for S2 (alcohol sensor), (a) sensitivity curve, (b) behaviour to 
temperature, and (c) step response 


Table 1. Specifications of the MOS sensors in the e-nose 


Sensor Hanwei Item Primary Gas Secondary Gases 
S1 MQ2 Combustible gas H2, LPG, CH, CO, Alcohol, Propane 
S2 MQ3 Alcohol CO, H2 
S3 MQ4 Methane Propane, Butane, Alcohol 
S4 MQ7 Carbon monoxide H2, LPG, CH, 
S5 MQ8 Hydrogen (ee) 
S6 MQ135 Air quality NH, NOx, Alcohol, Benzene, Smoke, CO2 
Controller ADCs 
S1}-— == 
S2 , 
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Figure 4. E-nose system: (a) block-level diagram and (b) electrical diagram of the air pump control 
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Figure 5. The activity diagram for the e-nose software 
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Figure 6. The sensors’ response during the acquisition process 


Once the 180 s acquisition process has concluded, the microcontroller gathers the sensors readings, 
converts the continuous data into digital form, and builds a data frame which is sent to the computer via 
USB. The computer further filters the samples with a ten-tap moving average filter to eliminate the 
undesirable peaks in the sensors’ responses. To obtain reliable plots with comparable maximum and 
minimum values, the sensors’ individual baselines were determined and eliminated using (1): 


Yi = Vi-Yio (1) 


where y; represents the output of the i"-sensor and yio its baseline. 

The final prototype is shown in Figure 7(a). Note that the microcontroller is fitted in the plastic base. 
The necessary connections between the sensor array and the electronic module are located inside. Figure 7(b) 
shows the implemented system. Inside the hermetic box, the e-nose (plastic base, sensor array, and electronic 
module) can be seen. Outside it, the air pump conveying the air samples can be appreciated. The air pump 
control circuitry and the power supply are also shown. The final prototype exhibits compact dimensions as 
shown in Figure 1(a); low mass (500 g), and low cost (150 USD). 


2.2. FDM analysis 

The filter diagonalization method (FDM) was initially developed for quantum dynamics calculations 
[27] and later used for nuclear magnetic resonance (NMR) signal processing [28], and leak detection in 
pipelines [29], [30]. It provides a nonlinear parametric method for time-domain signal analysis using a sum 
of damped sinusoids. FDM is traditionally used to solve harmonic inversion problems (HIP), delivering high- 
resolution spectra. Compared to Fourier analysis, it is not limited by the incertitude principle, thus providing 
high-quality spectra together with a high signal to noise ratio (SNR) without needing a large number of 
samples. 


Figure 7. First prototype is (a) control unit and sensors in the plastic base and (b) the prototype with the 
external air pump and the hermetic box 
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2.2.1. FDM formulation 

To explain the method, let us consider a complex unidimensional signal cy=c(nt) with values along 
with equidistant time intervals nt with n = 0, 1,..., N-1. The FDM aims to represent cn as a sum of damped 
sinusoids, as shown in (2), where wg = 27f;, — jy, are the complex frequencies of the signal, including the 
damping factor, and dx are the corresponding amplitudes. To solve (1), the FDM associates a correlation 
function described by the Hamiltonian operator 2, which has complex eigenvalues {wp}, thus cn can be 
further transformed into (3): 


K 
= —jntw 
Cc, = > d,e IN Ok 
k=1 


Cn = (Pole? by) 6) 


(2) 


The problem can be simplified to the diagonalization of the Hamiltonian operator Ê or, as discussed 
in [28], to the evolution operator 0 = e~/™. Briefly, a symmetric internal product operator defined by 
(a|b)=(bla) without the conjugate complex is used, where ®y is the initial state. Assuming that an orthonormal 
eigenvector set yg is used to perform diagonalization of the evolution operator as shown in (4): 


0 =P ug lhe) = Fe Fou YAY (4) 
k k 


and substituting (4) in (3), yields (5): 


dy = (Pole lWo) = (Y¥,|o)? (5) 


The resulting eigenvalues determine the position and width of the harmonics while the eigenvectors 
define their amplitudes and phases. Assuming a set created from the Krylov vectors, generated by the 


evolution operator ®, = U"@, = e~% ach ®, and according to (5), it yields: 

(®,|7 Pn) = (®,|Pm+1) = Cm+n+1 (6) 
as the set in non-orthonormal, the overlapping matrix can be calculated according to (7): 

(Py1Pn) = ("P| Po) = (P00 Go) = Canes (7) 


In (7) is strictly related to the values of the measured signal. Notation U’ can then be used, being 
this the overlapping matrix representation of dimensions M+1xM+1. Similarly, U! can be used for U. To 
reformulate (2), it is then necessary to solve the generalized eigenvalues problem as shown in (8): 


U'B, = u,U°B, (8) 
—jnwkT 


where ukg = e contains the lines of the spectrum and its corresponding widths. Eigenvectors Bk 
contain both amplitudes and phases. 


2.2.2. FDM implementation for the e-nose 
With the aim of analyzing the sensors’ experimental data, the FDM was implemented. The sensors’ 
readings were used as inputs Cn, and the FDM was used to estimate their spectra. The following steps were 
performed by the algorithm: 
— Taking into account the Nyquist criterion, the frequency interval [fmin fmax] in which the spectral analysis 
of signal c, will be performed is selected; c, is sampled at ff = Z 
— An angular frequency equidistant axis with values 22fmin<@j<2afnax, jJ=0, 1, 2,..., Kwin is created. Value 
Kwin is chosen as Kyin = MUmax fin) as suggested in the literature 


— Three symmetric complex matrices U® of dimensions Kwin X Kwin, with p=0, 1, 2 are determined. To 
calculate the elements that do not belong to the diagonal, (9) can be used, where f, and g, are the Fourier 
transforms of the first and second part of signal cn [31], detailed in (10). 
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ei? f(y’) — ei? f(y) + eg, (yp) +e? g,(0') 


(p) y= 
U (p, p ) ei? — eio (9) 
M 
O= e 
n=0 
2M 
gp) = Y a 
ae (10) 


In (11) is used to calculate the elements located in the diagonal: 


2M 
U® (9,9!) = Xu +1- |m- nee 
nzo (11) 


— Solve the generalized eigenvalues problem with (8), where the eigenvalues and eigenvectors are 
calculated using the QZ factorization algorithm [32]. 
— Select the complex amplitudes dx using (12): 


Kwin M 
qt’? pA Bye) CeIn; 
oe (12) 


— Use values œx and dx to estimate the spectrum with (13). Figure 8 shows a flowchart summarizing the 
FDM algorithm implementation. The following section will verify its performance and suitability for e- 
nose data processing. 
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Figure 8. FDM implementation algorithm for the e-nose 


3. RESULTS AND DISCUSSION 
3.1. E-nose data processing with FMD 

A series of experiments were designed using common gases to test the performance of the device, 
obtaining some preliminary results that are presented in this Section. The environmental conditions registered 
during the experiments were a temperature of 23 °C + 2 °C and relative humidity (RH) of 30%. For each 
experiment, the e-nose was subjected to a gas sample, and two features were obtained: i) the sensors response 
in the time domain consisted of a series of 180 samples for each of the six sensors; ii) the frequency spectrum 
for each sensor response was calculated using the FDM described in section 2.2. The spectrum of a sensor is 
given by (14): 


y= > annd(f — fen) (14) 
k 


where fi, and a; are the frequencies and amplitudes belonging to the k"-harmonic found in the response of the 
n"™-sensor, respectively. The system response to clean air was first verified. As expected, the sensor set 


Int J Elec & Comp Eng, Vol. 12, No. 2, April 2022: 1286-1298 


Int J Elec & Comp Eng ISSN: 2088-8708 O 1293 


showed constant low magnitude output values. The FDM-based spectrum revealed low amplitude peaks and 
frequencies with inexistent harmonics. It can be therefore concluded that the sensors’ response to clean air is 
negligible. Figure 9 shows the e-nose response to acetone (C3H6O). Note in Figure 9(a) how the sensors 
respond to the stimulus represented by the dotted line, especially sensors S4 and S2 followed by S6, this last 
indicating the air quality. In addition, a spectrum containing harmonics of relevant magnitude in a narrow 
frequency range can be observed as shown in Figure 9(b). Table 2 summarizes the frequency (Sn_F) and 
amplitude (Sn_A) values found by the FDM algorithm. 
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Figure 9. The e-nose response to acetone on (a) time response and (b) FDM-based spectrum 


Figures 10(a) and (b) show the time response and FDM-based spectrum to ethanol at 71.5%, 
respectively. As shown in Figure 10(a), sensors S2 and S4 exhibit the most significant responses; while S2 
shows the highest amplitude value, S4 shows the fastest rise time. Table 3 lists the harmonic values found by 
the FDM algorithm. Note the correspondence with Figure 10(b), values obtained from S2 and S4 outstand 
from the rest of the sensors. 
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Table 2. Acetone: amplitude (Sn_A) and frequency (Sn_F) values got for each sensor reading using FDM 


SILA SI F S2.A S2F SRA S3.F S4A S4F S5_A S5F S6 A _ SOF 
0.039 0.464 27.466 0.001 0.063 0.176 0.012 0433 2.859 0.003 4.580 0.014 
7.571 0.003 28.083 0.006 0.311 0.005 0.016 0410 1.024 0.015 1.061 0.030 
2.132 0.016 5.256 0.028 0.154 0.017 10.011 0 0.178 0.031 0.542 0.051 
0.545 0.037 1.195 0.045 0.025 0.030 4.978 0.012 0.166 0.048 0.584 0.063 
0.539 0.051 0.322 0.060 0.011 0.051 0.295 0.029 0.015 0.065 0.373 0.079 
0.041 0.072 0.110 0.082 0.761 0.042 0.027 0.082 0.206 0.090 
0.031 0.082 0.071 0.105 0.058 0.073 0.019 0.116 0.060 0.132 
0.011 0.121 0.041 0.122 0.020 0.087 0.016 0.151 0.201 0.153 

0.022 0.229 0.015 0.168 

0.015 0.197 

0.034 0.134 

0.014 0.172 

0.026 0.167 
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Figure 10. The e-nose response to ethanol on (a) time response and (b) FDM-based spectrum 
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Table 3. Ethanol: amplitude (Sn_A) and frequency (Sn_F) values got for each sensor reading using FOM 
SI A SIF S2.A S2F S3_A S3F S4A S4F S5A S5_F S6 A S6_F 
0.012 0.423 0.468 0.002 0.048 0.002 3.146 0.003 0.189 0.007 0.428 0.004 
1.852 0.005 0.393 0.011 0.025 0.012 1.487 0.017 0.109 0.023 0.248 0.017 


1.135 0.019 0.187 0.023 0.531 0.032 0.216 0.031 0.083 0.031 
0.421 0.032 0.071 0.037 0.489 0.047 0.029 0.042 0.024 0.048 
0.203 0.049 0.021 0.050 0.178 0.063 0.010 0.059 0.012 0.066 
0.071 0.067 0.071 0.079 

0.017 0.085 0.018 0.118 

0.108 0.182 0.032 0.129 

0.023 0.142 0.013 0.151 

0.015 0.161 


Finally, Figures 11(a) and (b) show the time response and FDM-based spectrum to gas butane 
(C4Hio), respectively. Note in Figure 11(a) that S4 exhibits the highest amplitude value. Table 4 shows the 
harmonics found by the FDM analysis. Note the correspondence with Figure 11(b), values obtained from S4 
outstand from the rest of the sensors. 


— == 56 
` === Odor Pulse 


Amplitude (V) 


0 20 40 60 80 100 120 140 160 180 
time (s) 


@) 


= 2 
0.05 0.1 0.15 0.2 0.25 


0 0.05 0.1 0.15 0.2 0.25 


ol 
R 
a 


0 0.05 0.1 d 0.15 0.2 


0 0.05 0.1 0.15 0.2 0.25 


0.05 0.1 0.15 0.2 0.25 


0.05 01 0.15 0.2 0.25 


frequency (Hz) 
(b) 


Figure 11. The e-nose response to gas butane on (a) time response and (b) FDM-based spectrum 
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Table 4. Gas butane: amplitude (Sn_A) and frequency (Sn_F) values got for each sensor reading using FDM 
SILA SI F S2.A S2F S3 A  S3_F S4A S4F S5A S5_F S6 A S6_F 
0.060 0.011 0.131 0.003 1.760 0.006 3.603 0.002 2.782 0.004 2.388 0.004 
0.467 0.008 0.098 0.015 1.136 0.017 1.191 0.013 2.049 0.019 1.322 0.017 
0.191 0.021 0.033 0.029 0.632 0.021 0.333 0.034 1.230 0.033 0.581 0.032 
0.017 0.046 0.012 0.197 0.175 0.045 0.056 0.057 0.677 0.046 0.496 0.046 
0.036 0.069 0.014 0.074 0.095 0.070 0.140 0.060 
0.015 0.115 0.016 0.240 
0.010 0.247 0.014 0.134 


3.2. Comparison with Fourier analysis 

The discrete Fourier transform (DFT) is the most commonly used method to solve the harmonic 
inversion problem (HIP) using the reliable and efficient fast Fourier transform (FFT) algorithm. However, the 
FFT is sensitive to time-frequency uncertainties which limit the resolution of the resulting spectrum. As the 
FFT spectrum resolution depends on the number of processed samples, the use of excessively large data sets 
is often necessary. In contrast, the FDM is a parametric method that, upon the use of linear algebra, extracts 
the parameters relevant for the construction of the signal spectrum. The FDM outperforms the FFT as it 
requires a lower amount of data to build the spectrum and does not restrain the spectrum resolution with 
uncertainties. 

Figure 12 compares the FDM and FFT spectra for the acetone sample. The comparison is limited to 
the signals of sensors S2 and S4, which show the clearest response to acetone. Note that the FDM produces 
more accurate high-resolution spectra with well-defined peaks. Harmonics are well separated among them, 
clearly showing the signal frequencies in the spectrum. In contrast, the FFT produces wide peaks that tend to 
merge between them; this low-resolution effect might lose or miss information on the spectrum harmonics. 
Note that the FDM-based spectrum exhibits a higher number of harmonics represented by well-defined 
peaks. It can be therefore concluded that the FDM-based spectrum is a clearer and more accurate 
representation for the e-nose signal processing. 

Current work focuses on using these spectra to perform gas classification. Approaches such as 
random forest [33] and convolutional neural networks (CNN) [34] are currently being considered. Future 
work will evaluate the possibility of optimizing the number of sensors while keeping satisfactory results and 
migrating the electronic unit to an field-programmable gate array (FPGA) system on chip (SoC) architecture 
[35]. Applications in wireless sensor networks (WSN) for domotics/inmotics [36] and automated farming 
[37] are foreseen. 
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Figure 12. Performance comparison between FDM and FFT (180 samples) 


4. CONCLUSION 

This research work has described a self-developed e-nose device implemented with low-cost 
components. It is devoted to detect potentially harmful volatile compounds in the environment. The prototype 
involves an array of six metal oxide semiconductor (MOS) sensors, an electronic module for data acquisition, 
and a computer-based information system for signal analysis and visualization. MOS sensors offer small size, 
low power consumption, fast response, and recovery times. Their response can be divided in three time- 
regions: the reference value (baseline), the rising time, and the resetting time. 

The hardware design and testing stages are now complete ensuring that the requirements have been 
fulfilled and that the prototype is functional. The filter diagonalization method (FDM) was implemented to 
calculate the harmonics involved in the acquired signals. The algorithm showed high precision using a low 
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number of samples. To our knowledge, FDM has not been previously explored for e-nose data processing 
and has the potential to become a valuable tool for this application. 
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